Avalanche dynamics and nonexponential relaxation 
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The theory of SOC, and related avalanche dynamics, is proposed as the origin of the ubiquitous 
nonexponential relaxation observed in complex systems. Introducing some scaling laws and relations 
we have obtained that the normalized relaxation function follows an stretched exponential decay and 
that the frequency spectrum follows an "1//" noise. Moreover, in the MF approach the relaxation 
is found to be exponential. 



I. INTRODUCTION 



Self-Organized criticality (SOC) was proposed as a 
general theory to understand the ubiquitous "1//" noise 
and self-similar fractal structures ||, in dynamical 
systems with many degrees of freedom. The SOC idea 
has been illustrated by a large variety of computer mod- 
els, including sandpiles |J, evolution jij], interface depin- 
ning H , and more ^ , in which a slow driving field leads 
to a stationary state with avalanches with widely dis- 
tributed sizes. However, most of the research has been 
devoted to the study of the stationary critical state. 

Another puzzle seeking a physical explanation is the 
empirical observation that relaxation phenomena in na- 
ture is not exponential. Nonexponential relaxation has 
been observed in different phenomena like mechanical, 
dielectric and magnetic relaxation, and in a large variety 
of materials, including amorphous polymers, supercooled 
liquids, glasses, spin glasses, and amorphous semiconduc- 
tors 0. These are complex systems where relaxation 
phenomena offer a large variety of challenging problems. 

Many theoretical models has been proposed to explain 
the ubiquitous nonexponential relaxation The 
physical origin of the nonexponential decay is attributed 
to the existence of fractal clusters M, cooperative dy- 
namics H, hierarchical constraints |lfj|] , or anomalous 
diffusion ]ll[|. All these models have their advantages 
and limitations and may be applied to some, but not all, 
experimental situations. 

In the present work we propose the SOC theory, and 
the related avalanche dynamics, as a physical mechanism 
leading to the nonexponential relaxation, often observed 
in nature. The present approach may be enclosed among 
those models which attribute the nonexponential behav- 
ior to the existence of fractal clusters, however the clus- 
ters (avalanches) in the SOC have an intrinsic dynamical 
nature, and correlations are present not only in space but 
also in time. 



II. CRITICAL RELAXATION 

In this section we analyze the relaxation of a self- 
organizing system near the critical state. We first start 
with a mean-field approach to the problem and later we 
analyze some scaling laws for the different magnitudes 
characterizing the dynamics of the system. 



A. Mean field theory 

The first step towards a comprehensive understanding 
of SOC and related phenomena is provided by mean-field 
(MF) theories, which give insight into the fundamental 
physical mechanism of the problem. The first MF analy- 
sis to sandpile models is due to Tang and Bak jl2|]. They 
mapped a critical height (Abelian) sandpile model in a 
two state cellular automaton and took the driving rate 
and the average height as control parameters. Later Cal- 
darelli, Maritan and Vendruscolo generalized the MF ap- 
proach of Tang and Bak in order to include stochastic 
rules Jl^] . On the other hand, Vespignani and Zapperi 
JLlj mapped the sandpile automata in a three state cel- 
lular automata on a d-dimensional lattice, and took the 
driving field h and the dissipation e rates as control pa- 
rameters. In this way, they obtained that the system is 
critical in the double limit h, e — > and h/e — > 0. The 
approach developed here is in the same line as that of 
Vespignani and Zapperi |l4[ ]. 

We divide the states that each site can assume in sta- 
ble, unstable and active, and denote their average densi- 
ties by p s , p u and p a , respectively. Stable sites are those 
that cannot become active by addition of grains. Unsta- 
ble sites are those that may become active by addition 
of grains. Active sites are relaxing and transfer energy 
to its nearest neighbors. Sites may receive energy either 
from the driving field or from its nearest active neigh- 
bors. The driving field is characterized by the driving 
rate h, the probability per unit time that a site will re- 
ceive a grain of energy from the driving field. On the 
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other hand, the probability per unit time that a site re- 
ceives a grain from nearest active neighbors is given by 
1 — (1 — pa) 9 — gp a , where g is the effective number of 
nearest neighbors and p a is small. Then, taking into ac- 
count both contributions, the local driving field is given 
by 



h\ oc = h 



9Pa 



(1) 



where p a °^ is the zero field average density of active 
sites and x the susceptibility, characterizing the linear 
response of the system to the external field. Substituting 
eq. (Q) in (||), comparing terms with similar order in h, 
and taking the physically admissible solutions, we obtain 
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This local field gives the probability per unit time that a 
site receive a grain, either form the driving field or from 
nearest active neighbors. 

Let u be the fraction of stable sites that become un- 
stable by addition of a grain, p the fraction of unstable 
sites that become active by addition of a grain, and q the 
fraction of active sites that become stable. With these 
definitions we write the following rate equations for the 
average densities 



where 



u(pg - 1) 



pq 



= 1 



(10) 



(11) 



Moreover, just at 9 = 9 C and h — > 0, from eq. (||) it 
follows that 



— p a =-[!-(!- q)ph loc ]p a + phlocPu , 



dt' 



qpa - uh\ oc p s , 



plus the normalization condition 



(2) 



(3) 



Pn 



-_h 1/2 . 



(12) 



Taking into account these results, from eqs. (||) and (Q) 
we obtain 



Ps = — + 0(h,\8 
ug 



(13) 



Ps + Pu + Pa = 1 



(4) 



The first term in the right hand side of eq. (g) charac- 
terizes the transition of active sites to stable or unstable 
sites. The correction (1 — q)ph\ oc takes into account that 
the fraction 1— q of active sites that become unstable may 
remain active by addition of grains from the local field. 
The second term of the same equation characterizes the 
transition of unstable sites to active induced by the local 
driving field. On the other hand, the first term in the 
right hand side of eq. (||) gives the fraction of active sites 
that become stable, and the second one, the transition of 
stable sites to unstable ones induced by the local driving 
field. 



After imposing stationarity 
eqs. (E[0) one obtains 



0, 



0), from 



upqgpl + [pq + u(l -pg+ pqh)] p a - uph = , (5) 



Pa 



uh + gp a 



Pu 



1 



Pa - Ps 



(6) 



(7) 



The first expression is a quadratic equation for p a . Its 
solution gives the average density of active sites as a func- 
tion of the driving field. We can expand p a {h) for small 
values of h, i.e. 



Pa= Pa 0) +Xh + O(h 2 ) 



(8) 
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(14) 



Thus, for 9 < 9 C , the zero field average density of ac- 
tive sites is zero while, for 9 > 9 C , it is proportional to 
(8 — C )P , with [3 = 1. The susceptibility diverges when 
9 — > 6> c , 9 7^ 9 C , according to \9 — # c | -7 , with 7 = 1. Just 
at 9 Cl the average density of active sites scales with the 
driving field as ft, 1 / 5 , with 5 — 2. These features are rem- 
iniscent of ordinary critical phenomena, h is the external 
field, 9 plays the role of temperature and p a is the order 
parameter. 

It is believed that conservation is a necessary condition 
to obtain SOC in sandpile models The global 

conservation law states that the average input flux hL d 
must balance, in average, the dissipated flux p a eL d . e is 
the dissipation rate per toppling event. It is an effective 
parameter that takes into account the dissipation, either 
in the bulk or at the boundary. Then, in the stationary 
state 



(15) 



From this equation one determines the average density 
of active sites as a function of the driving field, obtaining 
Pa = h/e. But, to be consistent with eq. (||), we must 

have = and x — l/ e - Moreover, using the expres- 
sion for the susceptibility obtained above, eq. (^?J), one 
obtains 



u 

6 . 

q 



(16) 
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In systems with dissipation in the bulk e is a fixed 
parameter and, therefore, the system is in a subcriti- 
cal regime 9 < 9 C . To obtain criticality (9 = 9 C ) we 
have to fine-tune e to zero. Hence, in this case, there is 
no difference with ordinary non-equilibrium critical phe- 
nomena, we have to fine-tune 9 (e) to reach the critical 
state. On the other hand, in systems with dissipation at 
the boundary, e decreases with increasing lattice site L. 
The number of boundary lattice sites, where dissipation 
may take place, growths slower than the total number of 
lattice sites. One can therefore assume that e ~ L~^, 
where /j is a scaling exponent. In the thermodynamic 
limit [L — ► oo) e = and, from eq. (jl|), it follows that 
9 = 9 C . In this case, the system self-organize itself into 
the critical state. Hence, the conservation law is not a 
sufficient condition to obtain SOC. The way in which 
energy is dissipated plays a determinant role. 

We have investigated the properties of the stationary 
state and the conditions to obtain criticality. Now we 
proceed to analyze the relaxation of the system after a 
switch off of the external field. If the system is close to 
equilibrium, h = and e ~ then form eqs. (|^), Jl4| ) 
and ( |l6| ) it follows that 



p(t) ~ exp(-pet). 



(17) 



Thus the relaxation towards equilibrium is exponential, 
with a relaxation time r e ~ e _1 . If the system is at the 
critical state (e = 0) then the system takes an infinite 
time to reach equilibrium. This phenomenon is char- 
acteristic of critical systems and it is known as critical 
relaxation. 



B. Scaling laws 

The MF theory leads to an exponential relaxation, 
however it is well known that relaxation in real systems 
is nonexponential. The spreading of an avalanche in MF 
theory can be described by a front consisting of non- 
interacting particles that can either trigger subsequent 
activity or die out. The MF theory thus neglects cor- 
relations among the particles. Next we investigate the 
propagation of a deltaic perturbation, in space and time, 
though a self-organizing system close to the critical state, 
but using some scaling laws. 

Let us calculate the average number of perturbed par- 
ticles tp(t) due to a deltaic perturbation, in space and 
time. ip(t) is then the response function of the system. 
These particles are causally connected in space and time, 
thus forming an avalanche with average size 



dtip(t) 



(18) 



We can express the average response if)(t) as a superpo- 
sition of the distribution of avalanche sizes, i.e. 



dsP(s)ip(s,t) , 



(19) 



where P(s) is the distribution of avalanche sizes and 
ij){s,t) is the reponse function of s-avalanches. In order 
to be consistent with eq. (18) ip(s,t) must satisfy 



dti/)(s, t) 



(20) 



Now we are going to assume scaling. If I is the lin- 
ear dimension of an s-avalanche then its characteristic 
size and duration should scale as l D and P, respectively, 
where D is the fractal dimension of the avalanches and 
z is a dynamic scaling exponent. Therefore, the charac- 
teristic time scales as s z l D and ip(s,t) may be written 
as 



V(s,t) = ^/i(fo-^) . 



(21) 



where f{x) is a cutoff function. The exponent q is de- 
termined from the normalization condition in eq. (pQ), 
obtaining q — D / z — 1. Moreover, tp(s, t) is independent 
of s for ts~ z l D < 1 then f(x) ~ etc for x <C 1. On the 
other hand, ip{s,t) -C 1 for ts~ z / D 1 and, therefore, 
f(x) < 1 for x > 1. 

The distribution of avalanche sizes near the SOC state 
satisfy the scaling relation 



(22) 



where r and a are scaling exponents and fa (x) is a cutoff 
function with similar asymptotic behavior as fi(x). e is 
the control parameter which determine how close is the 
system to the SOC state, i.e. the dissipation rate in the 
previous MF approach. 



Substituting eqs. ( |2l| ) and fl22[ ) in eq. ( |19| ) one obtains 
m = t^Mte 1 ^) , (23) 

where 



A = —a 

z 



D(2 - r) 



(24) 



and fs (x) is a cutoff function. 

Finally, the normalized relaxation function cf>(t) defined 
through the expression 



i t °°dt>m 



f dt>m ' 

may be approximated, for te 1 ^ <C 1, by 



exp 



(25) 



(26) 



where r e ~ e _1 / A . Thus, under the scaling laws assumed 
above the relaxation function follows the Kohlrausch 
stretched exponential. 
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The stretched exponent k is defined in eq. ([24]) through 
the scaling exponents r, D and z. However, using some 
scaling relations we can obtain a simpler expression. In 
systems at a SOC state the correlation length scales as 
the system size £ ~ L and £ ~ eT v ', where z/ is a scaling 
exponent. It is also possible to show that in homoge- 
neous SOC systems (s) ~ L 2 . Moreover, from eq. ( p2[ ) it 
follows that (s) - e-( 2 - T )/ '. Hence, 



(s) ~ L 2 ~ £ 2 ~ e~ 2v ~ e -( 2 " T )/ CT , 



(27) 



and therefore 



2 - t = 2ctza (28) 
On the other hand, the cutoff avalanche size is given by 



Sc ~ C 



and therefore 



LW = 1. 



(29) 



(30) 



Then, using (|24|), ( |28j ) and (|3C|), one obtains the scaling 
relations 



1 2 
A = — , k = - 

zv z 



(31) 



The stretched exponent k only depends on the dy- 
namic scaling exponent z (the duration of an avalanche 
scales with its linear dimension as P). In the MF ap- 
proach z = 2 in order to recover the exponential relax- 
ation obtained in the preceding subsection. However, if 
one considers correlations among the different branches 
forming an avalanche then one expects that the duration 
of the avalanche grows faster with its linear dimension 
and, therefore, z > 2. 



III. DISCUSSION 

We have obtained that the normalized relaxation func- 
tion follows an stretched exponential decay, with a 
stretched exponent given by the second scaling relation in 
eq. (|3l|). However, this functional dependence is limited 
to the time scale 1 < t < r e ~ e 1/A . In the SOC state 
e = and, therefore, the stretched exponential decay will 
be observed for infinitely long times. 

This temporal behavior leads to the frequency spec- 
trum 



S(f) 



dt cos(2?r/ 't)(f>(t) ~ /" 



(32) 



in the frequency range e -1 / A « / < 1. We thus found 
that the avalanche dynamics leads to a "1//" noise in 
the intermediate frequency window. This result seems to 
be similar to that reported by Bak, Tang and Wiesenfeld 
in their pioneering work M , however this is not the case. 



Bak, Tang and Wiesenfeld obtained a " 1/ /" noise assum- 
ing a power law distribution P(T) ~ f avalanche 
durations, which is equivalent to assume a power law dis- 
tribution of avalanches sizes, but they also assume that 
each avalanche relaxes exponentially, i.e. 



P(T) exp(t/T) - t 



-(a-l) 



(33) 



thus obtaining S(f) ~ / 2+a . The exponent a may be 
expressed in terms of of the scaling exponents introduced 
here. Since T ~ s z ' D and P(s)ds = P(T)dT then 

^, (34) 



a = 1 



D, 

- T-l =1 

Z 



where the second equality is obtained using the scaling 
relations in equations (|8|) and (|30|). 

In both cases the power spectrum is of the form 1//^, 
but (3 = 1 + K is larger than one in our model, while 
j3 = 1 — (D — 2)/z is smaller than one in the Bak, Tang 
and Wiesenfeld approach. Moreover, the functional de- 
pendence of the normalized relaxation function is differ- 
ent, we obtain a stretched exponential while they obtain 
a power law decay. This discrepancy is a consequence 
of the assumption made by Bak, Tang and Wiesenfeld 
|l| that the avalanches relax exponentially, with a relax- 
ation time equal to its duration. Such an assumption 
does not take into account the intrinsic dynamic nature 
of an avalanche, which may lead by itself to a nonexpo- 
nential relaxation. 

The normalized relaxation function associated with 
an avalanche of duration T, or equivalently of size s 
(T - s z /°), can be obtained from eqs. @) and (|2£ 
resulting 



4>(T,t) = 



g /T) dx X ^' D h{x) 



f~dxxi-'/ D fi(x) 



(35) 



Hence, the characteristic time is the avalanche duration, 
but the normalized relaxation function is not necessarily 
exponential. It is well known that cooperative dynamics 
1 eads, without assuming any disorder, to a nonexponen- 
tial relaxation dynamics. 



IV. CONCLUSIONS 

We have proposed the avalanche dynamics as the origin 
of the ubiquitous nonexponential relaxation observed in 
complex systems. Introducing some scaling laws and re- 
lations we have obtained that the normalized relaxation 
function follows an stretched exponential decay with an 
stretched exponent k = 2/z, where z is the dynamic scal- 
ing exponent. Moreover, in the MF approach n = 1 and 
the relaxation is exponential. 

The stretched exponential decay leads to the "1/ /" 
noise spectrum S(f) ~ / _1_K - This result is more gen- 
eral than the one reported by Bak, Tang, and Wiesenfeld 
0, since it takes into account the intrinsic dynamic na- 
ture of the avalaches. 
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